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Abstract 

Motivated by recent progress in experimental techniques of electric dipole mo- 
ment (EDM) measurements, we study correlations between the neutron and elec- 
tron EDMs in common super symmetric models. These include minimal supergrav- 
ity (mSUGRA) with small CP phases, mSUGRA with a heavy SUSY spectrum, 
the decoupling scenario and split SUSY. In most cases, the electron and neutron 
EDMs are found to be observable in the next round of EDM experiments. They 
exhibit certain correlation patterns. For example, if dn ^ 10^^^ e cm is found, de 
is predicted to lie in the range 10^^^ — 10~^^ e cm. 



1 Introduction 



The current EDM limits [1] - [3]^ 

\dn\ < 6 X 10~^^ ecm , 
141 < 2 X 10"^''' e cm , 

I^Hgl < 2 X 10~^^ ecm (1) 

impose severe constraints on CP violating physics beyond the Standard Model and, 
in particular, super symmetry. Super symmetric models accommodating the above 
bounds usually predict EDMs not far from the current limits which makes the next 
round of EDM experiments with the sensitivity d„ ~ 10~^^ e cm, dg ~ 10~^° e cm 
[4] - [8] particularly interesting. Other EDM experiments such as measurements 
of the deuteron and muon EDMs [9] would provide important complementary 
information on CP violation in supersymmetry. 

In this work, we address the questions "What are the EDM expectations in 
SUSY models ?" and "How can one distinguish different sources of EDMs ?" . It is 
well known that generic SUSY models predict too large EDMs which constitutes 
the SUSY CP problem [10]. The problem is resolved in certain classes of supersym- 
mctric models, of which we choose four representative types. These include SUSY 
models with 

• small CP phases 

• heavy spectrum 

• decoupling 

• split SUSY 

To determine the source of EDMs, we study correlations between the neutron 
and electron electric dipole moments. In particular, if the neutron EDM is due to 
the QCD ^-term [11] - [13], one expects very small leptonic EDMs. On the other 
hand, in supersymmetry both the liadronic and leptonic EDMs are enhanced and 
there exist certain correlation patterns. Determination of such patterns is the main 
subject of this work. 

The relevant low energy Lagrangian describing interactions of electrons, quarks, 
gluons and photons is given by 

£ = e^GG + f-bcG-&G' - Ldf 7(Fa)75/ - '^gsd, q{GaHq , 

where w, df, dq are the Weinberg operator [14] coefficient, the fermion EDM, and 
the quark chromo-EDM (CEDM), respectively. Here F and G represent the photon 
and gluon field strengths, respectively. The electron EDM is given by de, whereas 

^The EDM collaboration has recently announced an improved (preliminary) neutron EDM bound 
<3x 10-26 e cm [4]. 
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the neutron EDM is a model-dependent function of 9, w, dq and dq. We will use 
the naive dimensional analysis (NDA) approach [15],^ 

~ 2dd - O.bdu + e{OAdd - O.ldu) + 0.3 GeV x ew , (2) 

where we have set 6 = 0. Here the Wilson coefficients dq,dq,w are evaluated at 
the electroweak scale. The chromo-EDM and the Weinberg operator contributions 
involve considerable uncertainties. Furthermore, the result is sensitive to the quark 
masses which we choose as mu{Mz) = 2 MeV and md{Mz) = 4 MeV. The approach 
based on QCD sum rules gives a somewhat similar result [17], 

d^^~2dd-0.5(i„ + e(Jrf + 0.5d«) + 0.1 GeVxe'w . (3) 

In the Weinberg operator contribution, we have used the (model-dependent) renor- 
malization factor of Ref. [15]. In practice, the two approaches usually agree within 
a factor of 2 (unless w dominates) which suffices for our purposes. We note that 
there are also neutron models which include the strange quark contribution [18]. 
This effect is difficult to estimate and involves large uncertainties. We defer a study 
of such models until a subsequent publication. 

2 EDMs in the Standard Model 

We start by considering EDMs induced by the QCD ^-term. This is the primary 
source of hadronic EDMs in the Standard Model [11] - [13], [19]. In particular [13], 



d„ ~ 3 X 10"^^ 6* e cm (4) 

with about 50% uncertainty. On the other hand, the electroweak contributions are 
very small, d„ 10^'^^ e cm [20]. To satisfy the experimental bound (^P), the 6 
parameter has to be tiny, < 0(10^^*^). Such a small value can hardly be explained 
by the Standard Model, which constitutes the "strong CP problem". The most 
popular solution to this problem in extensions of the SM invokes an anomalous 
Peccei-Quinn symmetry [21], which sets 9 to zero. However, this symmetry is 
expected to be broken by higher dimensional operators generated at the Planck 
scale [22], such that the resulting 9 is finite but small. There are also alternative 
solutions to the strong CP problem which employ other symmetries [23]. In these 
models, a finite 6 can be induced by radiative corrections. In either case, a small 
neutron EDM is expected but it is difficult to make any quantitative prediction. 

The 0-term also induces nuclear EDMs. In particular, it generates EDMs of the 
deuteron and the mercury atom. The former receives contributions from both the 

^The QCD correction factor for the du,d contributions has recently been recalculated [16] and found 
to be a factor of 2 or so smaller than that in [15]. This leads to somewhat smaller estimates for dn- 
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constituent nucleons and nucleon interactions, and can be evaluated via QCD sum 
rules with about 50% uncertainty [24] . An EDM of the mercury atom is induced by 
the Schiff moment which appears due to CP violating isoscalar and isovector pion- 
nucleon couplings. Despite recent progress in evaluating these contributions [25], 
there is still an order of magnitude uncertainty in d}ig{9), whereas dug induced by 
the quark (colour ) EDM contributions is understood much better. The result is 

dz) ~ -1 X 10"^^ 6* e cm , 

I^Hgl ~ 0(10"^^ - 10"^^) e e cm. (5) 

The above formulae together with suppressed leptonic EDMs provide a correlation 
pattern for the 0— induced electric dipole moments. 

It is important to remember that de is not measured directly, but instead is 
derived from atomic EDMs. In particular, the current EEDM bound is due to 
the thallium atom EDM measurement [2] and the relation dxi — —585 de- If the 
6-teim is non-zero, this relation is altered. Although no reliable calculation of 
dTi{6) is available, rough estimates [26] of the valence proton EDM contribution 
give dxi ~ 10~^° 9 e cm. On the other hand, the SM electroweak interactions 
(including Majorana neutrinos) usually induce de of order lO^^^e cm [27]. Thus, in 
the Standard Model with 9 > 10~^^, dxi is dominated by the ^-background and is 
less sensitive to de- This illustrates that the effects of the 0-term can be important 
for atomic systems and should be taken into account. 

3 EDMs in supersymmetry 

In supersymmetric models, there are additional sources of CP violation associated 
with complex phases in the SUSY breaking F-tcrms and, in addition, flavour mis- 
alignment between the Yukawa matrices and the soft breaking terms [28]. EDMs 
are generated already at the one loop level and typically exceed the experimental 
bounds by orders of magnitude [10]. This constitutes the SUSY CP problem. The 
problem is alleviated in certain classes of supersymmetric models. These include 
models with small SUSY CP phases, those with a heavy SUSY spectrum, the 
decoupling scenario and split SUSY. In what follows, we study neutron-electron 
EDM correlations in these types of models. 

We note that there are also other possibilities for EDM suppression. For in- 
stance, CP violation may have flavour off-diagonal nature due to some symme- 
try [29], [30]. We defer a study of this option until a subsequent publication. Fi- 
nally, EDMs may be suppressed due to accidental cancellations among independent 
terms [31]. This option is however disfavoured by the mercury EDM constraint [32] 
- [35]. 

In our numerical analysis, we choose two representative values of tan (3, 5 and 
35, and analyze the NEDM-EEDM correlations separately for SUSY CP violation 

4 



Figure 1: One loop SUSY contributions to fermion (C)EDMs. 



induced by the phase of the /i-term, (p^, and CP violation due to the phase of the 
trilinear A-terms, (pA- We scan over the parameter space of a given model and 
present our results as scatter plots de vs dn- 

3.1 Small CP phases 

First, we study the minimal SUGRA model (mSUGRA) with smah (lO^^) Qp 
phases. The smallness of the physical phases may be due to approximate align- 
ment between the phases of the soft terms as occurs in the dilaton-domination 
scenario [36] with the Giudice-Masiero mechanism for the ^-term, or due to some 
approximate symmetry [37]. 

The model is defined in terms of the following GUT scale parameters, 

"T-o , "1-1/2 , 1^1 , tan/3 (6) 

and CP-phases cpA, (pfi- Here is the universal scalar mass, mi/2 is the universal 
gaugino mass and A is the trilinear parameter. The Higgs potential parameters 
I// 1 and are found by imposing radiative electroweak symmetry breaking. We 
assume that the CP phases associated with the gaugino mass and the S//-term 
have been rotated away by appropriate f/(l)R and C/(1)pq transformations. In this 
case, the physical CP phases are parametrized by (pA and (p^. 

The (C)EDMs are dominated by one-loop diagrams involving gluinos, charginos 
and neutralinos, Fig^ 

d, = dl + df+df , 
dq — c/g ~\~ ~\~ d^ , 

de = df + df . (7) 

The two loop contributions are considerably smaller. The relevant formulae can be 
found, for example, in Ibrahim and Nath, Ref. [31]. 

To get a feeling for the size of the EDMs, let us consider a simple approximation 
0^ ~ (/"A = 1, tan/3 ~ 3 and assume a single mass scale M for the SUSY 



parameters (at the electroweak scale). Then, 

/300 GeV\^ , _24 

I M I ^^^^^'^^ e cm , 

/300 GeV\^ , 25 

I M J sm<?i X lO^^^ecm . (8) 

Both dn and (ig grow hnearly with tan/3. Clearly, for the SUSY spectrum with 
electroweak masses, the CP phase has to be of order 10"^. We see that the neutron 
and electron EDMs differ by about an order of magnitude, however no prediction 
of their magnitudes can be made since the result depends on how small the phase 
is. If no de at the level 10"^'^ e cm is found, the CP phase will have to be smaller 
than 10~^ which appears highly unnatural. It would be rather difficult to engineer 
a robust mechanism which would force the phases in the Lagrangian to align with 
such an extraordinary precision. Thus, one may argue that de and dn in this scenario 
should not be far below the current experimental limits. 

In the same simple approximation, the mercury atom EDM is given by dug ~ 
(300 GeV/Mf sin^x IQ-^e g cm. Here dng is dominated by the quark CEDMs [32] 
and we have included an extra factor ^1/4 due to a recent reevaluation of the 
nuclear /atomic matrix elements. The deuteron EDM is expected to be similar to 
the neutron EDM, do ~ dn- 

Let now turn to our numerical results presented in Figs. 16171 In these plots, we 
vary mo,mi/2, \A\ randomly in the range 200 GeV - 1 TeV and the phases cpA, 
in the range [— 7r/500, vr/SOO]. In the left plot, (j)A is set to zero and, in the right plot, 
(j)^ = 0. For non-zero (j)^, there is a clear linear dg-dn correlation and de is about 
an order of magnitude below dn, as expected, de is dominated by the chargino 
diagram, whereas dn receives comparable contributions from the charginos and 
gluinos. At higher tan /?, the EDMs increase linearly. For 7^ 0, the correlation 
is less pronounced. For a given dn, the spread of de values is about two orders of 
magnitude. This is because d„ is dominated by the gluino diagrams, whereas dg 
is dominated by the neutralino contributions. The former are relatively insensitive 
to tuq since the squark masses are dominated by the gluino RG contributions, 
whereas the latter are sensitive to mo through both the slepton masses and the 
/i-parameter. Thus, fixing de does not determine d„ accurately. Note that EDMs 
induced by cpA do not receive tan/3-enhancement. 



3.2 Heavy SUSY spectrum 

The SUSY contributions to EDMs are suppressed if the entire SUSY spectrum is 
in the TeV range (Eq.®). Such a possibility is motivated by the strong bound on 
the Higgs mass which requires the stop mass to be of order 1 TeV. If all SUSY 
masses are controlled by the same scale, the spectrum is heavy. This senario can 
be motivated in various ways, see e.g. Ref. [38]. 
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Figure 2: 2 loop EDM contributions. The diagram on the left induces the Weinberg 
operator. On the right, is an example of the Barr-Zee type diagram (a is a pseudoscalar 
Higgs). 

In this class of models, the EDMs are usually dominated by the one loop di- 
agrams of Fig^ and the analysis is very similar to that presented in the previous 
subsection. We note that if the de experiments with the 10"^*^ e cm sensitivity yield 
a null result, it would imply that the scale of SUSY masses is 100 TeV (Eq.®). 
Models with such a high SUSY breaking scale are disfavoured by the gauge coupling 
unification and naturalness considerations. Thus, again one expects a non-zero re- 
sult in the next round of EDM experiments. 

In our numerical analysis, we study the mSUGRA model (Eq.©) with 
mo,mi/2, \A\ in the range 2 TeV - 10 TeV and (pA, i^i the range [— 7r,7r]. Our 
results are shown in Figs. 18191 The d^-de correlation patterns are very similar to 
those in the small phase scenario, namely there is a well defined correlation in the 
(f)^ ^ case, whereas for cpA 7^ it is far less pronounced. 

3.3 Decoupling 

The largest contributions to EDMs come from sfermions of the first two genera- 
tions. If these are very heavy, > 10 TeV, the most dangerous contributions are 
suppressed [39]. The third generation is required to be light by naturalness and 
contributes to EDMs at the 2 loop level. This decoupling scenario can be realized, 
for instance, in certain types of GUT models where the hierarchy between the first 
two and third generation masses appears due to RG running [40]. 

The neutron EDM is dominated by the stop and sbottom contributions to the 
Weinberg operator [41], whereas the electron EDM is due to the Barr-Zee [42] 
type 2 loop diagrams [43], Fig|5J This class of diagrams also includes graphs with 
internal charginos and h,H,A Higgs bosons [44]. Thus, 

dn ^ dn{w) , 

4 = ^Barr-Zee(/,x+,X°) . (9) 

We note that the Barr-Zee type diagrams also contribute to the neutron EDM, 
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Figure 3: (left) and de (right) in the decoupling scenario, (/)^ e ['I'/IO, 7r/2], tan/3 = 

3. The GUT scale parameters A — mi/2 — 



are in GeV and the EDMs are in e cm. 



The QCD sum rules model gives d^ a factor of a few smaller. 



but these are suppressed compared to the Weinberg operator [34]. The correlation 
between dg and d„ is rather subtle since (i„ depends on the gluino and the third 
generation squark masses, whereas de depends on the latter as well as the chargino 
and Higgs masses. An order of magnitude estimate of the resulting EDMs can be 
obtained by setting all SUSY masses (apart from the first two generation sfermions) 
to be M ^ Mz and the CP phases to be given by a single quantity (p. Then, for 
moderate tan/3 [41] - [44], 

2 



dr. 



dp. 



300 GeV 
M ) 

300 GeVV 



tan/3 sincj) x 10 '^^e cm 



M 



tan P sin ( 



10~^^e cm 



(10) 



Here dn and de differ by almost two orders of magnitude, whereas in other scenarios 
this hierarchy is usually a factor of ten. 

In this approximation, the mercury EDM is less sensitive to the CP phase (p. 
This is because the Weinberg operator contribution to dng is negligible and the 
EDM is generated by the subleading quark CEDMs. When M or (/> in Ea. ((Tn)) is ad- 
justed to satisfy the current neutron EDM bound, the resulting dug is < 0(10^^^) e 
cm, which is likely to be below the reach of the mercury EDM experiments. The 
deuteron EDM, however, is again found to be similar to dn- 

In our numerical analysis, we use the following GUT input parameters 



m 



1/2 ) 



1^1 , tan/3 



ill) 



and CP phases (pA, <Pfj.- Here is the mass parameter for all the scalars apart from 
those for the first two sfermion generations, which are assumed to be decoupled. 
As in the mSUGRA case, |//| and \B^\ are determined by radiative electroweak 
symmetry breaking. We vary m^,mii2, \A\ in the range 200 GeV - 1 TeV and 0^, 
(p^ in the range [— vr, vr]. Fig|31shows typical ranges of the dn and de values for order 



one CP phases . Clearly, d j2 10 e cm and d, 



10 



-27 



10 e cm are expected 



^dn and de are dominated by contributions sensitive to (p^. 
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in this scenario. Smaller EDMs would imply that either the CP phases are small 
or that the SUSY spectrum is heavy, including the third generation sfermions. 

The corresponding correlations are presented in Figs. llUlllI For (p^ 7^ 0, and 
dn can differ by two or one order of magnitude, depending on the balance between 
the Weinberg operator and the chargino Barr-Zee contributions. At larger tan/3, 
the situation remains the same except the EDMs increase proportionally. Note that 
since we scan over CP phases in the region [— 7r,7r], we include the possibility that 
the CP phases are small such that smaller values of the EDMs compared to those in 
FigElare allowed. For 7^ 0, the hierarchy between dn and de increases to three or 
four orders of magnitude. This is because dn is dominated by the stop contribution 
to the Weinberg operator, whereas the leading chargino Barr-Zee contribution to 
de is now absent and de is due to the stop Barr-Zee diagram. At larger tan/3, 
this hierarchy reduces since unlike the Weinberg operator, Barr-Zee contributions 
receive tan/3-enhancement. 



3.4 Split SUSY 

This is an extreme version of the decoupling scenario in which the third generation 
sfermions are decoupled as well and naturalness is abandoned [45], [46]. From the 
EDM perspective, it provides an interesting framework in which the neutron and 
electron EDMs are generated by the same type of diagrams and thus are highly 
correlated. 

The EDMs are induced by a version of the Barr-Zee type diagram, FigEl with 
the sfermion loop replaced by the chargino loop [44]. Since all the Higgses except 
for the SM-like h are assumed to be heavy, the chargino loop is attached to the 
fermion line through h and a photon [46], or and W~ [47].^ The relevant 
formulae can be found in Refs. [44], [47]. Clearly, the contributing diagrams are 
the same for the electron and the neutron such that the EEDM and NEDM are 
strongly correlated. We have therefore, 

^ _ jBarr-Zce(x+,X°) 

de = dBarr-Zee{x+,X°) . (12) 

Simple estimates are obtained by setting the chargino mass scale M to be much 
larger than the Higgs mass ~ 100 GeV [46], 

300 GeVy sin( 



•SOOGeVysin^ 

' M y tan^ ^ ' 



^It has recently been shown [48] that similar diagrams mediated by h and Z are also significant for 
the neutron EDM. This does not affect our numerical estimates. 
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M2 

Figure 4: d„ in split SUSY. Here (f)^, 4>m^ are varied in the range [-tt, tt]; e [200 

GeV, 1 TeV] and m/^ e [100 GeV, 300 GeV]. 

Here dn is obtained by rescaling with the factor mq/nie ~ 10. The tan (3 sup- 
pression can be traced down to the fact that the EDMs are due to the chargino 
mass matrix rephasing invariant M11AI22M12M21 oc sin/3 cos /3 and the SM-hke 
Higgs couphngs have no tan/3-enhancement. 

We note that since no CEDMs are generated at this level, the mercury EDM 
is suppressed. The deuteron EDM is, as usual, of order d„. 

In our numerical analysis, we vary the GUT scale gaugino masses Mi 2,3 in 
the range 200 GeV - 1 TeV and cpfj, in the range [— vr, vr]. We set the EW scale /i- 
parameter by hand in the range 200 GeV - 1 TeV and also vary ruh in the range 100 
GeV - 300 GeV. Typical dn values are shown in FigllJ whereas the corresponding 
de is found by a simple rescaling. For order one CP phases and low tan /3, dn is 
between 10~^^ and 10"^"^ e cm. Smaller values are obtained at large tan/3 or for 
small CP phases/heavy gauginos. 

The dn-de correlations are presented in FigEl As expected, dn is almost in 
one-to-one correspondence with d^,. The small spread in the EDM values comes 
from an interplay between the chargino and neutralino Barr-Zee diagrams, which 
contribute in different proportions to de and dn- In the left figure, the broad tail 
of dn below 10~^^ e cm is an artifact of our numerical procedure and appears due 
to 1-loop contributions of heavy but not completely decoupled sfermions (we took 
mn^mA = 10^ GeV and rrisferm = 10^ GeV). 

4 Conclusions 

We have analyzed correlations between the electron and the neutron EDMs in 
typical supersymmetric models. Unlike in the case of the 0-term induced electric 
dipole moments, in SUSY leptonic EDMs are enhanced and are likely to be observed 
in the next round of experiments. 

Assuming that all SUSY CP phases are order one at the GUT scale, lower 
bounds on the EDMs can be obtained in the decoupling and the spht SUSY sce- 



10 



narios: 



decoupling : 4 ~ (10^^ - 10"^) dn > 10"^^ e cm , 
split SUSY : de ~ 10"^ dn > 10-29 e cm , 

where the lower bound is saturated for a TeV range SUSY spectrum (apart from 
the first two generation sfermions) in the decoupling case, and at large tan /? in the 
case of split SUSY. 

For mSUGRA with small CP phases or a heavy spectrum, no solid lower bound 
can be derived. However, non-observation of de at the level 10"^'' e cm would imply 
that either the CP phases are < 10"^ or the scale of SUSY masses is > 100 TeV. 
Both of these options make super symmetric models very unappealing, so if low 
energy supersymmetry is indeed realized in nature, one expects dg > 10^'^'^ e cm. 

Concerning the dn-de correlations, the main feature of supersymmetry is that 
it enhances the leptonic EDMs and we observe the following correlations: 

small phases : de 
heavy spectrum : de 
decoupling : de 
split SUSY : de 

where we have assumed that all SUSY CP phases are of the same order of magni- 
tude. These relations as well as correlations with d^ and dng can help distinguish 
supersymmetry from other new physics models. 

We note that these correlations are quite stable to variations in the GUT scale 
soft breaking parameters. As long as de and d„ are generated by similar sets of 
diagrams, the most important factor for de/dn is me/mq and other effects being 
subleading.^ 

In this paper, we have focused on the low and moderate tan /? regimes where 
the effects of four-fermion operators (Fig.©) are less important [49] - [51] and 
have neglected certain two loop RG effects [52]. Also, we have not imposed other 
experimental constraints such as the bound on the Higgs mass, etc. which are 
expected to restrict the parameter space further. Clearly, such effects will not 
change the qualitative picture. 

Finally, we note that in the decoupling and the small phases scenarios, at least 
part of the SUSY spectrum is light and can be observed at the LHC. On the 
other hand, if the CP phases are order one, the entire SUSY spectrum may lie in 
the multi-TeV range. Then direct discovery of the superpartners at the LHC may 
prove to be quite challenging, whereas the EDMs can still probe such a possibility. 

^In particular, breaking the universality between the GUT scale squark and slepton masses does not 
affect our results significantly: the squark masses at low energies are dictated by the gluino mass and 
are much less sensitive to m^^^^j.^{GlJT) . 
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~ 10-1 dn , 

~ 10-1 dn , 

~ (IQ-l - 10-2) dn , 

~ 10-1 dn , 



e 



A 



H 



D 



D 



(a) 



(b) 



Figure 5: An example of a four-fermion operator inducing EDMs of composite objects 
[50]. Here, a pseudoscalar (a) and scalar (b) Higgs exchange induces a CP odd contact 
interaction dd ei-y^e contributing to rfxi- 

Acknowledgements. We are grateful to P. Harris and M. Pospelov for helpful 
communications. We would also like to thank D. Doyle for collaboration at the early 
stage of this project. 



12 





10-'' 




10-^' 




io-'« 






-a' 


10-^' 




10-" 




10-" 




e 10- 



d„ (ecm) 



10-'° 10-^' 10-'* 10-2' 10-'' 10-' 
d„ (ecm) 



Figure 6: de vs c?„ in mSUGRA with small phases, tan/? = 5. Left: 0^ 0, right: 
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Figure 7: As in FiglHI but for tan/3 = 35. 
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Figure 8: de vs dn in mSUGRA with a heavy spectrum, tan/3 = 5. Left: (/)p 7^ 0, right: 
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Figure 10: de vs d„ in the decoupling scenario, tan/3 = 5. Left: (f>^ ^ 0, right: 0a ^ 0. 
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Figure 11: As in FigHni but for tan/3 = 35. 
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Figure 12: de vs d„ in split SUSY. Left: tan/3 = 5, right: tan/3 = 35. 
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